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ABSTRACT 

QSOs from SDSS, 2QZ and 2SLAQ covering an order of magnitude in luminosity 
at fixed redshift exhibit similar amplitudes of clustering, with the brightest sample 
showing a clustering length only 11 ±9% higher than the faintest sample. In addition, 
QSO clustering evolution at z > 0.5 is well fitted by a model that assumes a fixed host 
halo mass. If halo and black-hole masses are related, then this may imply that QSOs 
occur in a relatively narrow range of halo masses with a correspondingly narrow range 
of BH mass. HST and Gemini high resolution imaging of QSOs covering a large range 
in luminosity also show a relatively narrow range in QSO host galaxy luminosity. We 
argue that the slow evolution of early-type galaxies out to z ~ 1 — 2 may also provide 
further support for a slow evolution of QSO host BH masses. The result would mean 
that if high- z QSOs radiate at Eddington rates then low-z Syl must radiate at « 100 x 
less than Eddington. We discuss the consequences in terms of four empirical models 
where (a) QSOs radiate at a fixed fraction of Lsdd, (b) QSO luminosity 'flickers' over 
time, (c) QSOs have a single BH mass and (d) QSOs are long-lived and evolve via 
Pure Luminosity Evolution (PLE). We conclude that the L>Edd model requires Mbh 
and Mhaio to be decoupled to circumvent the clustering results. While the single BH 
mass and flickering models fit the z > 0.5 clustering results, they appear to be rejected 
by the Mbh — L relation found from reverberation mapping at z « 0. We find that 
the inclusion of z < 0.5 QSO clustering data improves the fit of a long-lived QSO 
model and suggest that the predictions of the PLE model for QSO BH masses agree 
reasonably with UV-bump and reverberation estimates. 
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1 INTRODUCTION 

Over recent years there has been much new information on 
QSO clustering at high and low redshift. Here we combine 
results from various surveys to look specifically for the lu- 
minosity dependence of QSO clustering. 

The QSO clustering amplitude in the range 0.5 < z < 
2.2 evolves only slowly with redshift. Since the mass cluster- 
ing is falling as redshift increases range, this means that 
the QSO bias appears to rise with redshift. Converting 
bias into redshift via Press-Schechter theory implies that, 
for example, 2QZ QSO s at most red shifts inhabit haloes 
of > 2 x 10 12 M^ (|Croom et all 120051 ). By comparing the 
space density of QSOs with the space density of such haloes 
jMartini fc Weinberg 200 ll ). the general conclusion is that 
QSOs are short lived, with lifetimes of ~ 10 7 yrs. This short 
lifetime has seemed a necessary ingredient in a picture where 
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the QSO halo mass does not app ear to increase with time 
under gravity l|Croom et al.l [20051 ) . 

QSO halo mass and BH mass are generally as- 
su med to be correlat ed usin g the empirical re lation s 
of iFerrarese fc Merrittl d2000l); rGebhardt et"aH (|2000h : 
iFerraresei (|2002T l; IWvithe fc Loebl |2005l ) and "wTshall ini- 
tially also be making such an assumption. Som e empirical 
suppo rt is found in the observational analysis of iFine et al.l 
(2006) where BH masses estimated from BLR line widths 
are compared to halo masses estima ted from QSO clustering. 
The analyses of lCroom et alj (2005) suggest that it depends 
on which of the above halo-BH mass relations are assumed 
as to whether there is evolution in the Eddington ratio with 
redshift or not. Certainly the errors are large enough to en- 
compass various possibilities. 

ICroom et all (120051) found little evidence for luminos- 
ity dependent QSO clustering at fixed z in the 2QZ survey, 
however the luminosity range in any i ndividual magnitude 
limite d surv ey is necessarily lim ited. So lPorciani fc Norberd 
( 2006) and IShen et al l (|2009l ) claimed some evidence for 
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luminosity dependent QSO clustering respectively within 
the individual 2QZ and SDSS QSO surv eys but both re- 
sults w ere only marginally significant. Then lda Angela et al.l 
(2008) found little evidence of a luminosity dependence of 
clustering after adding the fainter 2SLAQ survey to the 2QZ 
survey. They suggested that t he results might be in good 
agreeme nt with the models of iLidz et all <|2006f > where the 
results of lCroom et al. I <|2005l ) were interpreted as suggesting 
that if QSOs have a short lifetime and evolve quickly from 
bright to faint states ('flickering'), then it might be expected 
that the luminosit y dependence of QS O clustering might be 
small as found bv lCroom et al.l |2005h . 

Here we shall first find further evidence for the in- 
dependence of QSO clustering and luminosity supporting 
these previous conclusions. We shall then consider the im- 
plications for QSO evolution by considering four simple 
empirical models where (a) Q SOs radiate at a fixed frac- 
tion Of LEdd f ee iReed Il984h. ( b) QSO luminosity 'flick- 
ers' over time l|Lidz et all 120061). (c) QSOs have a single 
BH mass l|Croom et al.ll2005l : Ida Angela et al.ll2008l ) and (d) 
QSOs are long-lived and follow Pure Luminosity Evolutio n 
(PLE, iBovle et al.lll987l ; iBovle. Shanks, fc Petersonlll988h . 
We shall attempt to test each of these against the cluster- 
ing and other QSO data. In the cases where models fail the 
QSO clustering test, we shall also consider whether the as- 
sumptions that underpin the clustering test may be at fault 
rather than the model. 

The data that we will use comes from the SDSS DR5, 
2QZ and 2SLAQ QSO redshift surveys. These surveys pro- 
vide a range of magnitude limits which cover a factor of 10 
in QSO luminosity. This means we can for the first time ef- 
ficiently check the dependence of small-scale clustering on 
QSO luminosity. We shall use these results as a context for 
a discussion on the phenomenology of QSO formation and 
evolution. 



2 QSO CLUSTERING DATA 

2.1 SDSS, 2QZ and 2SLAQ surveys 

Previously ICroom et all ([20051 ) used the 2QZ survey to es- 
timate the QSO correlation function and its dependence on 
redshift and luminosity. This was a survey of « 22655 QSOs 
in two w 375deg 2 strips in the NGC and SGC. The magni- 
tude limit was 18.25 < bj < 20.85 and the resulting QSO sky 
density was 31deg -2 . At the average redshift o f z = 1.4, the 
averag e absolute magnitude is Mj , ~ — 25mag. lCroom et al.l 
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measured s = bAt^Hh^Mpc and 7 = 1.2 ± 0.1 in 
the range 1 < s < 25h _1 Mpc for the amplitude and slope of 
the redshift-space co r relati on function, £(s). 

Ida Angela et al.l (120081 ) then used the 2SLAQ survey of 
9418 QSOs based on SDSS imaging to test the luminos- 
ity dependence of the QSO clustering. The magnitude limit 
was 20.85 < <7ab(~ bj) < 21.85 with a sky density of the 
resulting QSO sky density was w 49deg -2 , including the 
2dF QSOs where the two surveys overlapped, in a total area 
of 192deg 2 . The a verage absolute magnit ude at z = 1.4 is 
Mb j ~ -23.5mag. Ida Angela et al.l (120081 ) found a £(s) am- 
plitude and slope similar to that for 2QZ and again found 
that this QSO clustering amplitude depended very little on 
redshift or luminosity. 



Most recently. iRoss et al.l (|2009l ) have analysed the clus- 
tering of 30239 QSOs in the 5740deg 2 SDSS DR5 survey to 
a limit of Iab = 19.1, producing an average absolute mag- 
nitude of Mi « -26.3mag or M bj sss -26.1 atzxi 1.4. This 
provi des a sky density of 5.3deg~ 2 in the uniform sample 
where IRoss et all (|2009l ) measured s = 5.95 ± 0.45/i _1 Mpc 
and 7 = 1.16±o;ig in the 1 < s < 25h _1 Mpc range. 



2.2 Clustering comparison 

In Fig. 1 we now compare the clustering results from the 
three surveys directly with each other in terms of the 
redshift-space correlation function, £(s). The cosmology as- 
sumed in all cases is Qa = 0.7, fi m = 0.3. We see that all 
three results are quite consistent, particularly in the range 
5-20h _1 Mpc where the signal is excellent. The power-law 
amplitudes and slopes the authors have fitted are also con- 
sistent. We checked this by fitting the £(s) results consis- 
tently in the 1 < s < 30h _1 Mpc range. We have fitted for 
the real space correlation function scale- length, ro, assuming 
power-law slope, 7 = 1.8, infall parameter, ft = 0.4, and line- 
of-sight pairwise velocity dispersion < w 2 > 1 ' 2 = 750kms _1 
The latter, somewhat degenerate, parameters are in the 
range usually found in redshift-space distortion analys es of 
QSO surveys (|Hovle et al.ll2002l ; Ida Angela et al.ll200Sl ). We 
find good consistency in these results with SDSS giving 
r = 6.30±0.3h _1 Mpc, 2QZ giving r = 5.75±0.25h _1 Mpc 
and 2SLAQ giving r = 5.70±0.35h _1 Mpc. The best overall 
fit is r = 5.90 ± 0.14h _1 Mpc. The fits are simple \ 2 nts , 
with no account taken of covariance between the correlation 
function points, since previous results have shown that the 
sparse QSO sampling leaves them approxima tely indepen- 
dent, at least at the bin siz es we have used l|Croom et al.l 
120051 ; Ida Angela et all 120081 ). The brightest QSO sample 
scale length (SDSS) is therefore only 11±9% larger than the 
faintest sample. The overall \ 2 goodness of fit of all 3 surveys 
to the r = 5.90 ± 0.18h _1 Mpc model is 46.2 on 36 degrees- 
of-freedom, acceptable at the 1.75<r level. Thus in surveys 
ranging from magnitude 19.1 to 21.85 or a factor of « 10 
in luminosity ( — 26.1 < M bj < —23.5), the QSO clustering 
amplitudes are remarkably similar. Since the three surveys 
have very similar redshift distributions, this means that the 
amplitude of clustering at fixed redshift (z « 1.4) is approx- 
imately independent of QSO luminosity. We further note 
that other authors have measured the cross-clustering of low 
redshift QSOs with LRGs and have fou nd a similar lack of 
depen dence on luminosity. For example, iMountrichas et al] 
( 2009) found that z « 0.7 LRGs cross-correlated with SDSS, 
2QZ and 2SLAQ QSOs with similar amplitudes. Thus both 
auto- and cross-correlation analyses support the remarkable 
idea that intrinsically bright and faint QSOs show similar 
clustering environments. 



3 HALO AND BH MASSES VIA CLUSTERING 
EVOLUTION 

We next consider the halo m odel fit to the b ias- 
redsh ift relations show n taken from lCroom et all (|2005h for 
2QZ, Ida Angela et all l|2008l) for the 2SLAQ survey and 
IRoss et"aT1 ( 20091 ) for the SDSS sample. These results have 
shown that the bias at each redshift is compatible assuming 
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Figure 1. The QSO correlation function compared for the 
SDSS, 2QZ and 2SLAQ surveys. These surveys span a factor of ~ 
10 in luminosity at fixed rcdshift (z fa 1.4) but the clustering am- 
plitudes are very similar. The model assumes £(r) = (r/rg) 



1.8 



with line-of-sight velocity dispersion < w > 



i/2 = 



750kms and 



infall parameter f){z = 1.4) 
tt m = 0.3 and Q A = 0.7. 



0.4. The assumed cosmology is 



linear fluctuations and Gaussian density distributions with 
the QSOs occupying a single halo mass of w 3 x 1O 12 M0 
at all redshifts. If we also assu me that there i s a host halo 
mass - BH mass relation (e.g iFerraresel 120021 ). this means 
that QSOs may contain the same BH-mass at all redshifts. A 
model with fi xed halo mas s of Mh aio = 3.0 ± 0.38 x 10 12 M o 
as found by ICroom et al.l (|2005t ) for 2QZ is compared to 
the evolution of £20 in Fig. [2^, and to the overall QSO bias- 
redshift relation ui[2]d. The normalisati on for the mas s mode l 
is also chosen to be the same as that of lCroom etHI l|2005h . 
A fit based on the SDSS, 2QZ and 2SLAQ QSO datasets 
gives a fitted halo mass of M ha io = 3.27 ± 0.41 x 10 12 M© 
with rms — 2.0 x 1O 12 M0, consistent with the result for 
2QZ. It can be seen that the single halo mass model is a very 
good fit to these data, although we note that this model does 
not represent the ev olution of an individual QSO. Following 
ICroom et alJ (|2005h . their equation (24) w hich assumes the 
unevolving relation, Mbh ~ M^^ from IFerraresel (I2002T ) 
then gives Mbh ~ 5 x 10 s Mq, again approximately inde- 
pendent of redshift. We note that the model dependence of 
this mass i s larg e with the different Mbh — Mdmh models of 
IFerraresel (|2002T ) and the evolving model of IWvithe fc Loebl 
1 20051 ) predicting masses b etween 10 s M© and 10 10 M Q at 
fixed z (ICroom et al.ll2005h . Nevertheless the relative Mbh 
change with z is much less model dependent and this is more 
important for our purposes here. 

These conclusions only apply if we base our results on 
the SDSS, 2QZ and 2SLAQ samples. Including the IRAS 
AGN, SDSS AGN and Keck points only increases the re- 
sult to Mhaio = 3.74 ± 0.57 x 10 12 M Q but the single halo 
mass model is now marginally rejected by the data at the 
2-3oievel (x 2 = 42.8 on 26 dof, P < 0.025), with the highest 
residuals coming from the SDSS AGN and Keck points (see 
Figs. [2} and there is evidence for an increase in halo and 



hence black-hole masses as we move to lower redshift. For 
example, the halo mass corresponding to the SDSS AGN 
sample atzfs 0.13 is M ha io = 9.2 ± 1.7 x 1O 12 M . It might 
be argued that low redshift Seyfert 2 galaxies may be a dif- 
ferent class from the Type I QSOs that dominate at higher 
redshift but given their lower luminosities, it might be ex- 
pected that their clustering amplitude represents a lower 
limit to that of Seyfert I's at low redshift. In a unified pic- 
ture these AGN might be expected also to be representative 
of Seyfert I's. One caveat is that the SDSS AGN sample 
may also contain a population of LINERS which may not 
be comparable to Seyfert I in their clustering properties. 
Later we shall also argue that their m 20 x higher space 
density means that only a fraction are expected to be ob- 
scured Seyfert I's. We also note that lHickox et aD l|201ll ) has 
suggested that obscured QSOs, with ro = 6.0 ± 0.6h _1 Mpc, 
may show a higher clustering amplitude than unobscured 
QSOs, with ro = 5.3 ± 0.6h _1 Mpc, in their QSO samples. 
We postpone further discussion of these two lower redshift 
clustering points until Sections 7 and 8. 



4 QSO HOST GALAXIES VIA DIRECT 
IMAGING 

We now look at the empirical evidence for the mass and 
lum inosity of QSO hosts from h i gh re solution direct imag- 
ing. ISchade. Bovle. fc Letawskvl ( 2000T ) used HST to image 
X-ray QSOs out to z ~ 0.5, where imaging decomposi- 
tion is most reliable. They combined these data with those 
from similar observ ation s of higher luminos ity QSOs from 
iBahcall et all l| 19971 ) and lBovce et all l|l998T ). Fig. 13c from 
these authors shows a flat distribution of host galaxy lumi- 
nosity with QSO luminosity. This is particularly the case 
considering that low luminosity hosts of high luminosity 
QSOs are particularly difficult to detect. These authors also 
note that when a further decomposition into a bulge and 
disk is performed, there may then be more dependence of 
host bulge luminosity on QSO luminosity. The errors on the 
bulge luminosity become increasingly large and may provide 
an explanation for the wider scatter th at is seen in this re - 
lation. There are also observ ations (e.g. iMerloni et al.ll2010h 
and model predictions (e.g. iLamastra et al.l 2010l ) suggest- 
ing that the bulg e-BH mass relation m ay not apply at high 
redshift (see also iMcLure et ai"1l2006t ). The potentially dif- 
ferent results for bulge and galaxy luminosities must also 
be borne in mind because luminosity or more accurately 
disc stellar mass may better correlated with DM halo mass 
than bulge mass (e.g Kormendv. Bender, fc Cornelll l201ll ; 



iKormendv fc Benderil201ll ) 



This work w as extended to higher redshift by 
ICroom et ail l|2004h using Gemini ALTAIR+NIRI observa- 
tions of 10 luminous QSOs out to z ~ 2. Only one host 
galaxy was detected and the upper limits on the others were 
compatible with simple passive evolution models of the host 
from the present day. Models where the hosts evolved at 
the same rate as the QSO luminosities were rejected at high 
significance. Again if we assume that the bulge-BH mass re- 
lation applies at high redshift then these results continue 
to imply a constant BH mass for high-z QSOs. Similar re- 
sults have also bee n reported by other authors. For example, 
iFlovd et alJ (|2004r ). in HST observations, also found only an 
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Figure 2. (a) f 2 versus rcdshift for the SDSS, 2QZ 
and 2SLAQ QSO surveys, for an assumed = 0.7, 

f^m = 0.3 cosmo logy. Also s hown are the results for 
the SDSS AGN of IWake et al.l ll2004h , the IRAS AGN of 
iGeorgantopoulos k. Shanks! dl994T) and the Keck QSO-LBG cross- 
correlatio n result of | A delberg er fc S teidcl (2005) with £20 as re- 
ported by Ida Angela et al.l d2008f ). The solid line represents the 
model for fixed halo mass of 3 X 1O 12 M . The das hed li ne rep- 
resents the evolution of the long-lived model of lFrvl which 
best fits the data with bg(z = 0) = 1.42 ± 0.03. (b) As above for 
the linear QSO bias, bq, now calculated from the real-space QSO 
correlation functions. 



f» 2 mag range of galaxy host luminosity for an ~ 4 — 5mag 
range of QSO luminosity, with little evidence for a correla- 
tion between the two, except possibly in their upper enve- 
lope (see their Fig. 7, upper panel). Thus the results from 
QSO clustering analyses and from QSO direct imaging ap- 
pear to support the suggestion that QSO host galaxies may 
have a luminosity or, by extrapolation, stellar mass that is 
independent of QSO luminosity. 



5 IMMEDIATE IMPLICATIONS FOR QSO 
MODELS 

Based on the SDSS, 2QZ and 2SLAQ samples, we have found 
that either the clustering method is not appropriate or that 
all QSOs dwell in similar sized haloes. This could further 
imply that the broad-lined QSO population contain BHs of 
similar mass, approximately independent of QSO luminos- 
ity. The clustering results are powerful with the capability to 
rule out simple models. For example, the results appear in- 
compatible with the simple model where QSOs all radiate at 
a fixed fraction of Lsdd, unless there is no relation between 
Mbh and Mhaio- Many other models that make this un- 
derlying assumption will also have this pro b lem w ith these 
data. For example, Fig. 10 of iBonoli et al.l l|2009l 'l predicts 
that at z = 1.5, Mb < -24 QSOs give r = 7h _1 Mpc and 
Mb < —22.5 QSOs give rn = 5.25h _1 Mpc, ie a difference of 
Ar = 1.75h _1 Mpc over this « 1.5 mag luminosity range. 
In comparison, over a wider 2.75 mag range, the data in Fig. 
[T] only gives Aro = 0.6 ± 0.46h _1 Mpc. Thus the predicted 
Aro = 1.75h _1 Mpc is marginally rejected at 2.5a. (We ig- 
nore the error on the model at fixed z because the predicted 
luminosity dependence clearly persists at all 2.) Similar con- 
clusions apply to the simplest v ersion of the semi-analytic 
model of iFanidakis et all |201ll ) for the evolution of black 
hole growth. Ot her empiric a l mod els such as the QSO 'flick- 
ering' model o f iLidz et al.l d2006h and the single BH mass 
QSO model of ICroom et all i|2005r ); Ida Angela et all (|2008h 
do better, since they were designed to explain the early in- 
dications of the luminosity independence of QSO clustering. 

Thus, QSO host halo and BH masses at fixed L/L* 
may be independent of redshift despite the fact that the 
QSO luminosity brightens by a factor of ~ 30 in the range 
< 2 < 2.2. Moreover, since QSOs at fixed redshift have the 
same clustering strength independent of luminosity or L/L* , 
this could imply that all QSOs have a narrow distribution 
of BH mass, independent of both redshift and luminosity, 
with QSOs becoming increasingly sub-Eddington at lower 
luminosities. 

This appears unlikely. However there are some pieces of 
evidence that this might be the case. First, the luminosity 
function of the early-type galaxies thought to host QSOs 
appear to evolve very little over the redshift range out to z ~ 
1-2 l|Metcalfe et al.ll2006l ; IWake et al.ll2006l : iBrown et all 
2007). If the LF and stellar mass function do not evolve 
much then it is a short step to supposing that the dynamics 
of the early-type do not evolve much either (although see 
iDaddi et aLlbOOSl ). Then it could be that the BH mass might 
also be unevolving. The biggest BH lie in the biggest bulges 
and this relation is then predicted not to evolve with time. 
The QSOs then have to inhabit a narrow range of bulge/BH 
mass but with a large range of luminosity at fixed 2. 

An escape from this conclusion could be that the high-2 
QSOs have moved along the galaxy luminosity /mass func- 
tion to higher galaxy luminosities/masses and hence have 
higher BH masses more consistent with being at Edding- 
ton. However, the fact that it has been so difficult to iden- 
tify QSO ho s ts at 2 ~ 2 in high resolution imaging led 
ICroom et al.l f2004) to suggest that the data are consistent 
with the hosts being similar luminosity to those at lower 
redshift. 

Another possible argument for a slowly changing BH 
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mass with time is the fact that the velocity widths, 
of the broad lines are similar for high luminos- 



< v > 



1/2 



ity QSOs at high z as they are for low-luminosity Seyfert 
Fs at low z. Now in the standard view where QSOs ra- 
diate at approximately a fixed fraction of the Eddington 
luminosity, this is just a coincidence where the increase 
in BH mass by a factor of ~ 30 at high z is approx- 
imately offset by a comparable increase in the BLR ra- 
dius. Since the BL R radius-luminosity relation takes the 
form r sa L°- 5 -°- 7 dPeterson et al.l 12004 iKaspi et alj|2005l ; 
iBentz et al.ll2009T ) , if this relation is unevolving with redshift 
then < v 2 >oc L/r cx L 0,5-0 3 and so < v 2 >^ 2 oc L 25 " 015 , 
perhaps slow enough to explain the similarity in line widths. 
In this fixed fraction of LEdd case, there is therefore some 
empirical evidence to help explain this coincidence. Never- 
theless, the explanation is still vulnerable if the r — L relation 
actually proves to evolve with redshift. 

Alternatively, to maintain the single BH mass model, 
we could pursue the idea motivated by the early-type galax- 
ies that nothing evolves except the QSO luminosity. On this 
view, the BLR radius might then remain constant with lu- 
minosity and hence redshift ie r ~ L°. A fixed broad- line 
velocity width, < v 2 > 1//2 , and radius, r, gives a fixed BH 
mass with luminosity /redshift. The question then is whether 
it is feasible, for example, that the factor of w 30 increase in 
L* over < z < 2.2 could leave the BLR radius unchanged 
if a balance has to be maintained between photoionisation 
and gravity and this seems physically unlikely. It also seems 
at odds with the above evidence for the r — L (and Mbh — L) 
correlation found from reverberation mapping. Thus if these 
empirical results are correct then the constancy of QSO ve- 
locity width with L and z must remain simply a coincidence, 
rather than providing further evidence for the single BH 
mass model. 

Indeed, any single BH mass/flickering model motivated 
by QSO clustering clearly has a severe problem with the 
QSO Mbh — L relation fo und from reverberatio n mapping 
at low redshift. Equn (9) of lPeterson et al. (2004) rejects no 
correlation ie a slope of zero in the Mbh — L relation at the 
8.8a level! So here we have potentially the most serious and 
direct conflict with any simple single black hole mass inter- 
pretat ion of the QSO clustering results. As iPeterson et al.l 
l|2004h point out, the simple fixed fraction of Eddington 
model fits the reverberation data very well but this model 
most clearly contradicts the clustering results. If we are to 
reconcile the clustering data with the reverberation data, we 
must conclude that there must be some decoupling of Mbh 
and Mdmh and this immediately takes away the evidence 
for the flickering and any other single BH mass model. Mean- 
while, we postpone our discussion of reverberation mapping 
and the PLE model to Section 8 and Fig. U where we shall 
directly compare BH masses predicted by the PLE model to 
masses derived from reverberation mapping. 



6 NUMBER DENSITIES OF QSOS AT ALL 
REDSHIFTS AND LUMINOSITIES. 

The space number density of QSOs compared to the num- 
bers of remnant nuclear black holes in early-type galaxies 
has b een used to arg u e for short (r tti 10 7 yr) QSO life- 
times [Richstone^etjiL fl998). However, if the QSOs were 



restricted to a small range of host galaxy/halo/BH masses as 
implied by the similarity of the QSO clustering results then 
this argument would start to push QSO lifetimes higher. 
Generally, these estimates depend on an assumption that 
QSO luminosity is proportional to host BH mass and this 
assumption is also challenge d by the lack of a lumin osity de- 
pendence of QSO clustering. iRichstone et alj (|l998l ) suggest 
that the bright QSO space density at z w 2 is ~ 0.1% of 
the galaxy density at z — 0. Thus, for example, if the QSOs 
are restricted to « 1% of galaxies then this lifetime would 
increase to « 10 9 yrs. Clearly the more restricted the galaxy 
luminosity range, the longer these duty cycle arguments al- 
low the QSO lifetime to be. 

More detailed arguments for short QSO lifetimes have 
been similarly made on the basis of high pe a ks mo dels of 
biased QSO clustering. iMartini fe Weinberg! (|200lh made 
predictions for the QSO correlation function based on the 
halo number densities inferred from QSO lifetimes estimated 
from QSO and model halo densities. Using the ACDM re- 
sults from th e se aut hors and the 2QZ clustering amplitude, 
ICroom et all < |2005f > derived a QSO lifetime of 4-50 Myr at 
z ~ 2. We note that a monotonic increase of halo mass with 
QSO luminosity is an implicit assumption here, otherwise it 
is impossible to establish a minimum halo mass at a fixed 
QSO luminosity limit. In what follows we shall also make 
this assumption. If we take the change in b as m 9% from 
2SLAQ (r = 5.7h -1 Mpc) to SDSS QSOs (r = G.Sh^Upc) 
then fro m the approximate re l ation, b cx M^ 2 n , fitted to 
Fig. 6 of IMartini fc Weinberg! (|200ll ). this corresponds to 
an ~ 1.5 x bigger halo mass for the ACDM model. Now 
assuming an approximate ACDM halo mass function with 
N(< M) cx M 015 for M < lO 14 h _1 M (|jenkins et al.ll200ll ) 
we find that the number of haloes in the 2 — 3 x 10 12 h _1 MQ 
mass range is « 10% of the total with M > 2 x lO 12 h _1 M . 
This means that the estimates of QSO lifetime would now 
lie in the range « 40 — 500 Myr. Since the clustering of the 
2SLAQ, 2QZ and SDSS samples are statistically consistent, 
the increase in QSO lifetime could even be larger than this 
estimate. 

This conclusion from the clustering is stronger than 
from the direct QSO imaging results. Here there is an 
« 2mag range in host luminosity or a factor of w 6, cor- 
responding to « 40% of the haloes with M > 2 X 10 12 M© 
giving an ~ 2 — 3x increase in QSO lifetime. Again, since the 
range for the host luminosities is clearly an upper limit due 
to no account being taken of galaxy magnitude errors, there 
is no necessary disagreement with the clustering results and 
indeed the QSO lifetimes could be much higher if as ex- 
pected the galaxy magnitude errors dominate the ~ 2mag 
scatter. We therefore suggest that the short QSO lifetimes 
previo usly derived on the basis of the IMartini fc Weinberg! 
i|200lT) arguments are at best lower limits if the QSOs are 
restricted to a range of halo masses and even on the basis of 
the present clustering data could imply that this lower limit 
is r > 10 8 yr. 



7 LONG-LIVED QSO CLUSTERING MODEL 

If QSO lifetimes become longer then there is a clear pre- 
dictio n for QS O clustering via th e long -lived model of iFrvl 
(|1996T ) (see also lCroom fc Shanks! (|l996h l. Here QSOs (and 
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their host galaxies) are assumed to form at fixed bias at high 
redshift and then only evolve as separate entities clustering 
under gravity. This gives the relation:- 

b Q (z) = l + [b Q (z = 0)-l]/D(z) 

where D(z) is the gravitational growth factor and D(z) — 
1/(1 + 2:) for the EdS model. This then produces the evolu- 
tion for £20 shown in Fig. The relation produces a sig- 
nificantly worse fit than the single halo mass model dis- 
cussed earlier, falling with redshift whereas the data appears 
to increase. If we now include the lowe r redshift re s ults o f 
iGeorgantopoulos fc Shanks! l|l994) and IWake et all (|2004r ) 
then the higher clustering implied at lower redshifts pro- 
duces improved agreement with the long-lived model predic- 
tion. We find that the best fit bias is b Q (z = 0) = 1.42±0.03. 
However, the model is still rejected at high significance with 
X 2 = 80.2 for 27 data points (26 d e grees of freedom). Note 
that following ICroom et al l (|2005h this bias for the long- 
lived model is uncorrected for non-linear effects. 

Given that the model assumes a constant bias for all 
QSOs at formation, this model intrinsically predicts a cor- 
relation function compatible with a constant halo mass for 
QSOs independent of luminosity, at least at fixed redshift. 
However, it is true that if there was a strong luminosity de- 
pendence at fixed z shown by the QSO data then this might 
also have been compatible with the long-lived model. 

We note that the cosmology also plays a role here. If 
we assume an EdS cosmology then £20 reduces by a factor 
~ 2 by z ~ 2. While the situation in respect of the single 
halo mass model remains the same as for the ACDM case, 
the situation in respect of the long-lived model improves 
significantly, with the model now fitting the high SDSS AGN 
clustering amplitude while being able to fit better the £20 
points at high redshift. The best fit bias is bg(z = 0) = 
1.71 ± 0.03, assuming a rCDM normalisation of £20 = 0.12. 
Overall the model is acceptable in a chi-square test at the 
« 8 % significance level (x 2 = 36.7 for 27 data points). 

ICroom et al.l (120051 ^) argued against a long-lived model 
on the grounds that accretion must take place and the halo 
masses would change contradicting the evidence from the 
constant mass clustering with redshift. However, this ar- 
gument is now weakened by the fact that if applied to 
early-type galaxies and LRGs a contradiction also results 
in the context of the ACDM cosmology which is avoided 
by appealing to halo occupation distributions (HODs). In- 
deed, the evolution of the luminosity function and cluster- 
ing of early-type galaxies themselves appear to be reason- 
ably compatible with long-lived PLE models. The LRG lu- 
minosity functio n appears consist ent with passive evolution 
out to 2 « 0.7 (|Wake et al.ll2006l V Indeed, the NIR H and 
K galaxy number counts, dominated by early-type galax- 
ies at K < 20, show little evidence of stro ng evolution out 
to z > 1 [Metcalfe et al.l l|2006l . l200ll . I1996T ). At large scales, 
IWake et al.l (|2008T l found that for separations r > lh x Mpc, 
the evolution of LRG clustering was reasonably consistent 
with the long-lived, biased evolution out to z ~ 0.6. At 
r < lh~ Mpc there was some evidence that the evolutionary 
model of 113 (l996) was too strong for the data. However, a 
HOD with a low galaxy merger ra te of only 2.4%/Gyr can fit 
the data. ISawangwit et all l|2009h also found that a simple 
virialised cluster model, which clearly also conserves galaxy 



numbers, can improve the fit of a long-lived model to the 
small-scale clustering at r < lh _1 Mpc. Although these re- 
sults generally only apply to 2 < 1, the fit of a HOD model 
with only a small amount of galaxy merging but with a halo 
merging rate which is c ompatible with ACD M suggests that 
the simple arguments of ICroom et al.l (2005) against similar 
long-lived models for QSOs may not be as strong as initially 
envisaged. 



8 PLE MODEL PREDICTIONS 

We next connect the above clustering results to the major 
observation of the evolution of the QSO luminosity func- 
tion. One model which predicts a 'long-lived' evolution with 
redshift is a pure luminosity evolution (PLE) model. If 
the dependence of QSO luminosity on redshift is given by 
L(z) = L(z — 0)(1 + z) 1 , then the dependence of BH mass 
on redshift is given by 



Mbh = 



L(z = 0) 



2(7-3/2) exc 2 



to(3.2 



.7-3/2 



(l + 2)^~ 3/2 ) 



where e is the BH efficiency. At the characteristic L* knee in 
the luminosity function, an M B (z = 0) = —21.5 — 51og h 70 
mag QSO has luminosity L*(z = 0) = 3.3 x 10 i4 h^ 2 ergs -1 
in the range 1 200- 10000 A (jMarshall et al.ll 1984 ) . We have as- 
sumed the EdS model for analytic simplicity but with fixed 
Universal age, io = 13.7 Gyr. For zero initial mass at 2 = 2.2, 
e = 0.4 and 7 = 3 the QSO then accretes Mbh « 10 9 M© by 
2 = 0. We have seen that the low redshift QSO clustering re- 
sults may prefer such a long-lived model. On the assumption 
that the halo and long-lived models are self-consistent, the 
rise in predicted clustering amplitudes from 2 ~ 2 to 2 ~ 
implies an increase of rss 10 x in halo mass, since b oc M®' 2 lo 
from Section 6. Thus if the halo and galaxy growth rates are 
not decoupled as assumed in the previous section, then Mbh 
would increase by a factor of « 100 x (since Mbh oc Ml^ 5 
from Section 3) and the PLE model would be rejected on 
the basis of its much slower predicted Mbh growth rate (see 
Fig. [3]). Even if the long-lived clustering model now gives 
an improved statistical fit with the inclusion of the low-z 
data in Figs. [2] the theoretical PLE Mbh growth rate and 
the theoretical long-lived model growth rate are mutually 
incompatible if Mbh and Mdmh follow the above relation. 

Fig. U shows the BH mass for 62 bright Palomar-Green 
QS Os as estim ated from fitting the 'UV bump' in their spec- 
tra (Laor 1990), compared to the BH masses estimated from 
the PLE equation above. At the high mass end, the PLE 
masse s seem in good agreement with the masses of lLaor] 
(1990), although at low BH masses, the PLE masses ap- 
pear too high. Overall the agreement seems reasonable but 
this is the result for a PLE model with negligible initial BH 
mass. A single BH mass model would depend on the scat- 
ter around ~ 1O 9 M0 and the errors on the UV bum p BH 
masse s to explain the range of BH mass shown by the lLaorl 
1 19901 ) estimates in Fig. H 

In a PLE model, the QSO is long-lived, thus explaining 
why the QSO space den s ity at L* remains the same with 
redshift jMarshalll Il985l; iBovle. Shanks, fc Peterson! Il988l ; 
ICroom et all I2OO4I 2009). In any other model this repre- 
sents a coincidence. Basically the same l uminosity functio n 
behaviour is also seen in X-ray surveys <|Bovle et alj|l993h . 
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- L(z=0)=3.3x10 44 ergs-' 
L(z=0)=3.3x10 45 ergs-' 



Figure 3. M BH (z) for PLE models for an L* QSO and a 10L* 
QSO, assuming 7 = 3, e = 0.4 and M BH (z = 2.2) = 1 X 1O 9 M . 
The L* model shows a relatively slow variation with z. The 10L* 
model shows a somewhat faster relation with z and a significant 
difference with the L* model at fixed z. 




UV bump 
Reverberation 



M BH (Reverb/UV bump)/10 M 

Figure 4. PLE masses from the Mbh(z) rel ation i n the text 
compared to the 'UV-bump' fitted masses of iLaod l|l990ll for 
62 Palomar -Green survey QSOs and reverberati on masses o f 35 
QSOs from lPeterson et al.l l|2004l l. The masses of lLaoJ l|l99dl ') as- 
sume /i75 = 1 and go = 0,5, s imilar to what is assumed for the 
masses of lPeterson et al.l 11200411 and for the PLE model. The PLE 
model also assumes 7 = 3, e = 0.4 and zero initial BH mass. A 
line of equal mass and a line where the PLE masses are 3x larger 
than the other masses are shown. 



There have been some reports that the X-ray LF is incon- 
sistent with PLE (jHasinger et alfeOQ j;lGiacconi et al.ll200l 
lAlexander et aljlioo^ ; IWorslev et al.ll20o4 ) but these devia- 
tions are generally small and at low Lx and in a broad-brush 
way at least, PLE remains an approximate fit to the X-ray 
LF. PLE has always been an excelle nt phenomenological fit 
to the optical QSO LF (although see iBongiorno et aLll2007l ) 



but the question has been whether it has any physical mean- 
ing. If exactly true the PLE model would suggest that low 
z Seyfert I's would have as high a BH mass as high-z QSOs 
at fixed L/L* because the BH fuelling rate at low redshift 
is small compared to that at hig h redshift (Marshalj |l985l ; 
iBovle, Shanks, fc Peterson]! 19881 ). 

A single BH mass independent of QSO luminosity is 
not, however, a prediction of PLE at fixed redshift. The 
equation above implies that for fixed 7 and e a QSO 10 x 
brighter than L* will produce a 10 x bigger black hole mass 
of 1O 1O M by z = 0. Indeed, for a 10L* QSO, from z = 2.2 
to z — 1.5, of order the minimum QSO lifetime, a black-hole 
mass of 5 x 10 9 Mq will already have been created. If we 
assume that the QSO has to be radiating at no more than 
Eddington, this implies that the initial mass for such QSOs 
has to be 1O 9 M0. If the fainter QSOs have a similar initial 
mass then this will minimise any Mbh difference between 
bright and faint QSOs (see Fig. O. 

Of course, if halo and BH masses are completely decou- 
pled, then there may be no need to accommodate the ap- 
proximate luminosity independence of clustering in this way. 
However, it could be that although halo/BH mass growth 
rates are decoupled, initial halo and BH masses could still 
be correlated. If so, the Mbh '■ Mhaio relation must evolve 
with redshift in the sense that the halo mass containing a 
particular Mbh must decrease as z increases. At z = 2.2, 
all QSOs would then start with the same initial Mbh and 
have the same halo mass to be consistent with the cluster- 
ing observations. Higher luminosity QSOs then accrete at 
a faster rate as in Fig. [3] but no luminosity dependence of 
the QSO clustering is generated because of the decoupling 
of halo and QSO/galaxy merging rates and the clustering of 
all QSOs grows at the long-lived rate. 

In the context of the original PLE model, QSOs at all 
luminosities at z ~ 2 could radiate at close to the Eddington 
limit and those at lower redshift then would radiate at in- 
creasingly sub-Eddington rates, up to « 30 x. In the above 
PLE model case where all the initial Mbh must now be 
the same to fit the clustering results, if we assume that at 
z ~ 2 the brightest QSO is radiating at Eddington, then at 
the same redshift the fainter QSOs have now to radiate at 
sa lOx lower than the Eddington rate. At z ~ 0, these intrin- 
sically faint QSOs have thus now to be radiating at « 300 x 
below Eddington. This may be physically feasible in that 
standard optically thick, physically thin, accretion discs can 
accrete at down to ~ 1% of the Eddington rate (C. Done, 
priv. comm.). Below this rate the accretion disc is unstable 
to evaporation. Thus the -27 mag QSOs at high redshift and 
the -22 mag Seyfert Ts locally can both be powered by thin 
disc accretion. Thus there appears to be a reasonable range 
of QSO luminosities where the PLE or singie-BH mass mod- 
els can operate within the context of a thin disc accretion 
model. 

The mean BH mass measu red in low luminosity SDSS 
AGN bv iHeckman et all (|2004l ) via the Mbh '■ o~ relation is 
w 1O 8 M0 and this seems in less good agree ment with th e 
» 10 9 M Q PLE prediction than the results of lLaon l|l990h . 
However, as we have already noted, the SDSS AGN com- 
prise some 18% of the galaxy population whereas broad 
lined QSOs comprise < 1%. Thus it is plausible that, al- 
though some fraction, (< 30%?) of the SDSS AGN may 
be obscured counterparts of QSOs, the remainder might be 
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expected to be a genuinely lower luminosity AGN popula- 
tion, sufficiently sub-Eddington to make any thin accretion 
disk unstable. Thus the obscured QSO population in a PLE 
model might lie only in the high BH mass tail of the SDSS 
AGN mass distribut ion which extends bey ond 1O 9 M0 (see 
Fig. 1 (left panel) of iHeckman et al l l|2004l ))- 

We also r evisit the reverberation BH masses of 
IPeterson et al.l (|2004T ) for 35 low redshift QSOs and com- 
pare these with their predic t ed PL E masses as we did for the 
'UV-bump' masses of lLaorl <|l99d l (see Fig. EJ. We see that 
there is a good correlation but that the reverberation masses 
are generally fa 5x higher than the PLE predic tions and th e 
agreement is less good than for the masses of lLaorl l|l990l) . 
For the 11 QSOs in common between the two samples, the 
'UV-bump' masses are as 3x higher than the reverberation 
masses but there is a significant scatter. The amplitude of 
the reverberation masses contains a scaling factor which was 
adjusted to fit the M bh ■ o"* relation in quiescent galaxies 
i|Peterson et al.ll2004) . We also note that if we chose to as- 
sume a fa 1O 9 M0 initial BH mass then the disagreement 
would clearly worsen between the PLE and reverberation 
masses. Nevertheless, the correlation with the reverberation 
masses seen in Fig. [4] means that the PLE model performs 
better in this respect than if the single BH mass or flickering 
models were assumed to apply at z as as at z > 0.5 (see 
Section 5). 



9 QSO X-RAY VARIABILITY 

A traditional problem for the PLE model is whether the 
massive BH it predicts, particularly in low redshift Seyfert 
Ps, can accommodate the short timescale variability ob- 
served in some of these sources. The timescale of variability 
set by a w 1O 9 M0 BH mass is « 10 4 s or ~ 3hrs. Now vari- 
ability in many QSOs is seen on timescales of < lhr and 
the question is whether these can be localised flares, affect- 
ing less than 5-10% of the area and total luminosity of the 
QSO. Clearly the smaller the fraction of the QSO luminos- 
ity that the flares affect, the more plausible the argument 
for the high PLE BH mass becomes. But the fact that the 
optical and X-ray evolutions are similar is a further coinci- 
dence that a long-lived model like PLE explains, since the 
optical accretion disc is expected to heat the X-ray 'corona'. 
We now consider these issues in more detail. 

The PLE and single BH mass models both make simple 
predictions for QSO variability. PLE predicts that at fixed 
L/L* , the amplitude of variability will be approximately in- 
dependent of redshift because the BH mass only changes 
slowly with redshift. The single BH mass model makes the 
even simpler prediction that the variability amplitude will 
be approximate ly independen t of re dshift and luminosity, 
i.e. L/L* or L. lAlmaini et al.l (|2000h found in a statistical 
X-ray variabi lity analysis of 86 QSOs f rom t he Deep ROSAT 
QSO survey iGeorgantopoulos et al.l l|l996h that generally 
the variability of QSOs shows a relatively flat distribution 
with redshift (see their Fig. 6b) . They also found that there 
was some evidence for an inverse correlation with Lx at 
z < 0.5 which was simi l ar to t hat found for low-z AGN by 
lLawrence fc Papadakisl (|l993l ); lNandra et all i|l997l ) but this 
correlation disappeared for z > 0.5 (see their Figs. 8a, b). 
More recently, other authors have made similar studies of 



QSOs in XMM and Chandra observations. In the CDF-S, 
iPaolillo et all (|2004T ) analysed 74 variable QSOs and found 
little evidence of dependence on Lx (see their Fig. 12). 
In the XMM - Newt on Deep Survey of the Lockman Hole, 
iMateos et"ai1 1)20071 ) found variances for 74 QSOs and AGN 
that were approximately indepen dent of luminosit y and r ed- 
shift (see their Fig. 5). Finally, IPapad akis et alJ (|2008h in 
the same XMM-Newton Lockman Hole QSO dataset found 
no significant correlations with either redshift or luminosity 
(see their Fig. 2). We conclude that while at z < 0.5, there 
seems evidence for an inverse correlation between luminosity 
and variability amplitude, at z > 0.5 there seems little evi- 
dence of any variability correlation with luminosity or red- 
shift. While the latter results argue for both the PLE and 
single BH mass model, the former result would only be con- 
sistent with the PLE hypothesis. Finally, we note that the 
amplitude of QSO variability is in the range a ~ 0.1 — 0.3. 
Thus this is in the range 10-30% of flux at which level it 
may be easier to argue that the high BH masses predicted 
by PLE models may be allowed, if only a fraction of the 
QSO luminosity is involved in the variability. 



10 DISCUSSION 

SDSS, 2QZ and 2SLAQ QSOs show similar clustering scale 
lengths to within ~ ±10%, while spanning an order of mag- 
nitude in luminosity. The QSOs in these surveys also show 
slow clustering evolution with redshift and a model with a 
single halo/BH mass for all QSOS fits both these results 
well. If true, then it implies that if the most luminous SDSS 
QSOs are radiating at Eddington luminosities at z ~ 2, 
then the faintest 2SLAQ QSOs at lower redshift must be 
radiating at an w 100 x lower rate than Eddington. Al- 
though problems with BH masses estimated from emission- 
line widths can be avoided by appealing to evolution in the 
BLR radius-luminosity relation or in the line-widths them- 
selves, it may be more difficult to explain away the strong 
Mbh — L correlation implied at z — from reverberation 
mapping. This presents an immediate problem for any single 
BH mass model and variants de signed to fit QSO clustering 
such as the 'flickering' model of iLidz et all (|2006l V 

If QSO halo mass is restricted to a very small range 
then estimates of QSO lifetime must rise significantly and 
our very rough estimate i s by a factor of ~ 40 t o ~ 1 9 yrs 
based on the arg uments of lMartini fc Weinberg! |200ll ). This 
then prompts considera tion of th e clustering prediction of a 
long-lived QSO model (|Fry||l996r ). We find that this model 
predicts an increase in clustering amplitude at z < 0.5. 
When we in c lude t he lo w-redshift AGN clustering results o f 
IWake et al.l (|2004h and IGeorgantopoulos fc Shanks! (|l994|) . 
we find improved agreement with the long-lived model. In- 
cluding these results also tends to worsen the fit of the single 
mass halo model. If the QSOs are long-lived then this also 
tends to argue against the idea that QSOs always radiate 
at a fixed fraction of LEdd because QSO BH mass tends 
to increase with time while QSO luminosity drops, so even 
if QSOs start at the Eddington rate they quickly become 
sub-Eddington. 

If the long-lived QSO model is correct, then the QSO 
clustering results may be consistent with the fit of PLE 
models to the QSO luminosity function. A physical inter- 
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pretation of why the QSO space density at L* is observed 
to remain approximately constant from z = 2 to z = is 
that the QSOs have a long lifetime and the luminous Type I 
QSOs seen at z = 2.2 dim by a factor of « 30 to become the 
low-luminosity Seyfert I's seen at z=0. Thus this PLE model 
prediction may be supported by the QSO clustering results. 
We have noted that if Mbh and Mdmh follow relations 
such as Mbh oc Md M hi then the theoretical growth rates 
of Mbh from PLE and the long-lived model are mutually 
incompatible. However, a key feature of the long-lived model 
is that it implicitly assumes that galaxy/BH and halo mass 
growth rates are decoupled so that the BH mass growth rate 
with 2 predicted by PLE is easily accommodated. At fixed 
z, PLE also predicts that BH mass growth rates are propor- 
tional to QSO luminosity. But again, since in the context of 
long-lived models, halo and BH mass growth are naturally 
decoupled there is no necessary contradiction with the clus- 
tering results here for PLE. This decoupling has given us 
the freedom to assume a constant initial BH mass for the 
PLE model so that it is actually in line with the luminosity 
independent clustering results at high z. Then the higher 
accretion rate of the more luminous QSOs gives the possi- 
bility of generating the luminosity-BH reverberation mass 
correlation found at z = 0, without also generating cluster- 
ing that is too luminosity dependent at these redshifts. In 
this picture, it is the lack of a relation between halo and 
galaxy/BH mass that explains the lack of luminosity depen- 
dence of QSO clustering seen in Fig. [1] rather than the lack 

of a relation bet ween luminosity and BH mass. 

Recently, IKormendy, Bender, fc Cornelll l|201ll ); 
IKormendy fc Bender! (|201 ll ) have presented evidence that 
nuclear velocity dispersions may not correlate well with 
galaxy circular velocities in cases they have observed where 
there is no bulge, just nuclear star-clusters. They also argue 
that dynamically measured Mbh do not correlate well with 
disc stellar masses. Since circular velocities and disk stellar 
masses usually correlate well with DM halo masses, they 
then conclude that DMH masses don't correlate well with 
Mbh- These results have some resonance with the argu- 
ments we have presented above based on QSO clustering 
and also the direct imaging of QSO hosts. We have argued 
that these results imply either that QSOs have the same 
BH mass or that the halo and BH mass growth rates are 
decoupled. In particular, our long-lived PLE model argues 
for a decoupling of halo and BH mass at least at z < 2.2 
to reconcile the slow evolution of QSO clustering with the 
fa ster predicted evo lution of halo clustering. In the models 
of lLidz et al.l (|2006l ) it is the QSO luminosity that decouples 
from the BH/DMH masses to explain the clustering results. 
But in these scenarios the BH and DMH masses are 
expected to be s trongly correlated, in contradi c tion with 
the new results of IKormendy. Bender. fc Cornelll (|201lT ) . 



11 CONCLUSIONS 

For clarity, we summarise the arguments for and against the 
four QSO models we have discussed above: 

• QSOs radiate at a fixed fraction of LEdd'- 
-in favour - reverberation mapping Mbh — L relation 
- agrees with virial estimates of Mbh using unevolved 
BLR L — r relation. 



- against - this model predicts strongly luminos- 
ity dependent clustering which is not observed; need 
to break the Mbh — Mdm h rela tion as suggested by 
IKormendy. Bender, fc Cornelll |201ll) . 

- if the QSOs are long-lived as implied by the luminos- 
ity independence of the clustering then they cannot always 
operate at a fixed fraction of LEdd and be consistent with 
observed QSO luminosity evolution. 

• The 'flickering' model of iLidz et all ^OOrj ): 

- in favour - L independent z > 0.5 QSO clustering. 

- against - the reverberation mapping Mbh — L relation 
at z ~ 0. 

- virial estimates of Mbh using the unevolved BLR L — r 
relation. 

- the longer duty cycles suggested by the coherent QSO 
LF evolution. 

- IKormendy. Bender, fc Cornelll i|201ll ) evidence against 
the assumed Mbh — Mdmh relation. 

• The single BH mass model: 

- in favour - L independent z > 0.5 QSO clustering. 

- against - the rise in clustering amplitude at z < 0.5. 

- the reverberation mapping Mbh — L relation at z ~ 0. 

- virial estimates of Mbh using the unevolved BLR radius- 
luminosity relation. 

• The long-lived/PLE model: 

- in favour - L independence of QSO clustering implies 
long lifetime. 

- not ruled out by the QSO clustering evolution over the 
< z < 2.2 range 

- reasonably reproduces reverberation and UV bump 
Mbh- 

- explains coherent QSO LF evolution. 

- invokes a decoupled Mbh — Mdmh r elation as evidenced 
by IKormendy. Bender, fc Cornelll (|201lf ). 

- against - virial estimates of Mbh using the unevolved 
BLR L — r relation. 

- short X-ray variability timescales for low-z Seyfert I. 

It is clear that models which predict well the QSO clus- 
tering results then tend to have a problem with the rever- 
beration and UV bump Mbh — L relation and vice versa. 
We take the view that it is easier to circumvent the clus- 
tering results (by decoupling Mbh and Mhalo) than the re- 
verberation mapping results. Therefore the answer to the 
question posed in the title may be that QSOs do not all 
have the same black-hole mass, although they may occupy 
similar halo masses at fixed redshift. Overall, we conclude 
that a model that is reasonably consistent with both cluster- 
ing and the reverberation/UV-bump masses is the long-lived 
PLE model. 
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